Transition radiation by matter-wave solitons in optical lattices 

A.V. Yulin, D.V. Skryabin, and P.St. J. Russell 

Department of Physics, University of Bath, Bath BA2 7AY, United Kingdom 

We demonstrate that matter-wave solitary pulses formed from Bose condensed atoms moving 
inside optical lattices continuously radiate dispersive matter waves with prescribed momentum. Our 
analytical results for the radiation parameters and the soliton decay rate are found to be in excellent 
agreement with numerical modelling performed for experimentally relevant parameters. 
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Recent observations of matter- wave solitons [1,2] have 
clearly been amongst the most breaking achievements in 
the burgeoning field of Bose-Einstein condensation of di- 
lute atomic gases. Balance between the spatial dispersion 
of matter waves and repulsive or attractive interatomic 
interactions ensures existence of dark [1] or bright [2] soli- 
tons, respectively. Dispersion of the atomic condensates 
can, however, be reversed by embedding the condensate 
into a periodic potential created by standing light waves, 
i.e. optical lattice [3-5]. The idea of changing the dis- 
persion sign and of the possible observation of the bright 
matter- wave solitons in the condensates with repulsive in- 
teratomic interaction has been around for a while, see [7] 
and references therein, and more theoretical results have 
been produced recently, see, e.g., [8-10], in the view of 
the rapid maturing of the experimental techniques [3-5] . 
The concept of the dispersion control by periodic poten- 
tials is also well known in solid-state physics [11] and a 
very active topic of research in nonlinear optics, see, e.g., 
[12]. 

To understand the initial motivation which leads to 
the results described below, it is instructive to recall 
the effect of the transition radiation known from clas- 
sical electrodynamics [13]. Transition radiation is a con- 
tinuous emission of electromagnetic waves by a charged 
particle moving with a constant velocity in spatially in- 
homogeneous medium. This radiation is emitted because 
the field created by the particle has different character- 
istics in different parts of the medium. When the par- 
ticle moves, the field reorganizes itself continuously and 
shakes off some of its parts in form of the radiation. Wc 
expect that similar phenomenon should take place with 
matter-wave solitons in optical lattices. 

In order to change the dispersion sign of matter waves 
forming the initially resting one-dimensional packet em- 
bedded into the optical lattice, one needs to position it in 
momentum space somewhere between the inflexion point 
of the energy-momentum, i.e., dispersion, characteristic 
and the edge of the Brillouin zone, see Fig. 1(a) [14]. 
Only exactly at the edge of the zone does the group 
velocity go to zero. Solitons with a spread of quasi- 
momenta centered at the edge will therefore be the only 
resting bright solitons in the condensates with repulsive 
inter-atomic interaction. Our primary interest below is, 
however, moving solitons. Once a solitonic wave packet 
moves through the periodic potential one can expect that 



its structure will not be able to instantly readjust itself 
to perfectly fit the conditions that the local density max- 
ima are positioned at the center of the local potential 
minima. Therefore the moving soliton has to continu- 
ously overcome the energy barrier, which leads to the 
continuous emission of matter waves. This energy barrier 
is analogous to the Peierls-Nabarro potential known for 
the solitons in discrete systems [6]. Using analogy with 
the electromagnetic case we term this effect as transition 
radiation of matter waves. 

We start our analysis from the Gross-Pitaevskii (GP) 
equation describing evolution of the macroscopic wave 
function of the zero-temperature Bose-Einstein conden- 
sate (BEC) interacting with off-resonant standing light 
wave. We assume that the condensate is tightly confined 
by the external harmonic potential along the Y and Z di- 
rections having the trap frequencies ujy.z = 2ir x 400s 
and that any deviations of the potential along X-axis 
from the sin 2 kiX produced by the intensity of the stand- 
ing laser field with wavenumber ki can be disregarded. 
We take for our estimates that ki = 27r/800n?7?.~ 1 and 
consider BEC made of 87 Rb atoms with a two-body 
scattering length a ~ 5.4 x 10~ 9 m. The character- 
istic transverse width of this BEC is then given by 
w = yf%J (mwy) ~ 0.5/im. Assuming that the conden- 
sate profile along the Y, Z directions is given by the low- 
est mode of the harmonic potential, we can derive the 
one-dimensional GP equation, which describes dynamics 
of the X dependent part of the full wave function. The 
dimcnsionlcss normalized form of this equation is [9] : 



id t yj — —d x ijj — Pyj cos 2a; + |V'| 2 V'- 



(1) 



Here the dimensionless time t and spatial coordinate 
x are measured respectively in the units of T — 
2m/(hkf) ~ 5 x 10~ 5 s and 1/fe; ~ 0.13/wn. Mean- 
ing of the dimensionless parameter (3 is easily inferred 
from the expression for the lattice potential in the phys- 
ical units, which is taken as 4(3E r sin 2 k[X , where E r — 
h 2 kf/(2m) ~ h x 20kHz is the recoil energy. To esti- 
mate number of atoms N in the condensate we introduce 
the effective area in the (Y, Z)-plane, A e ff = 2irw 2 ~ 
1.5/im 2 , and the atom density n ~ 10 14 cm -3 . Then 
N = J \rb\ 2 dx x A 2 eff km 2 / 3 /(8ir\a\) ~ 4 x 10 3 x / \yj\ 2 dx. 

We proceed by expanding -0 over the Bloch functions, 
b(x,k) [11]: ip(x,t) = J dkip(k,t)b(x,k). Hereb(x,k) are 
eigenfunctions of the operator L = d 2 + (3 cos 2x, such 
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that Lb = —e(k)b, where k is the quasi-momentum intro- 
duced as b(x,k) = g(x,k)e lxh [11], e(fc) is the energy of 
the non-interacting, i.e. linear, matter-waves, and g(x, k) 
is the function with the spatial period n. The first al- 
lowed, first forbidden (gap) and a small part of the 2nd 
allowed energy bands are shown in the energy-momentum 
plot in Fig. 1(a). 

Choosing k = k s e [—1,1], we expand g(x,k) in a 
Taylor series around k = k s and demonstrate that 

i>(x,t) = e ik ° x G s A(x,t). (2) 

Here, A(x,t) = J dk^(k,t)e^ k - k ^ x and G s 
is the linear differential operator: G s = 
^^ (l/n!)<9£(7(x, k s )(— id x ) n . Similarly Lip can be rep- 
resented as 

U = -e ik ° x G s £ s [-id x ]A(x, t), (3) 

where £,[-id x ] = £~ = o(l/m!)0j?e(fc a )(-i0 x ) ,n is the 
energy operator. Subscript V refers to quantities and 
functions calculated for k — k s . 

After substitution of Eqs. (2), (3) into Eq. (1), we then 
replace G s with its 1st order approximation g{x,k s ) = g s 
and £ s with £ 2s = c s — ie' s d x — l/2e"d x . Here e' s and 
e" are, respectively, group velocity and group velocity 
dispersion of matter-waves. Assuming smallness of the 
nonlincarity and averaging the resulting equation over 
g s we derive the renowned nonlinear Schrodinger (NLS) 
equation idtA = £2 S A + a s \A\ 2 A, where a s = J dx\g s \ . 
Below we are interested in the bright solitons with e s be- 
longing to the first allowed energy band. These are given 
by A s = R s (0 cxp{— int — ie s t}, where 
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R ^ = y— sech \N—\>z = x - £ ^ ( 4 ) 

sgnn — sgna s — — sgne" and n ^ is the nonlinearity 
induced energy shift. Details of the derivation of NLS 
equation from Eq. (1) has been previously published 
using the method of multiple scales [15]. Approxima- 
tions made above and in [15] imply that ip has sufficiently 
narrow spread of quasi-momenta around k = k s . How- 
ever, the approach introduced here is readily adaptable 
to give an access to the small amplitude corrections hav- 
ing quasi-momenta detuned far from k s , sec Eqs. (8), (9) 
below. Note, here that mobile envelope solitons (4) arc 
very different in their properties from practically immo- 
bile solitons occupying primarily one or few lattice sites 
and considered, e.g., in [9]. 

To introduce the effect studied and explained in this 
work, we first present the results of numerical mod- 
elling of Eq. (1) with initial conditions in the form 
g s {x)R s {x)e lksX for values of k s G (fco,l], e"(feo) = 0, 
ensuring that the effective mass, 1/e", is negative, see 
Fig. 1(a). Taking k s = 1, i.e. fixing soliton parame- 
ters at the point corresponding to the zero of the group 



velocity, we observe formation of the ideal solitary pulse 
[8] . The Fourier spectrum of this solution contains series 
of the equidistant peaks, the location of which is deter- 
mined by the spectrum of the corresponding Bloch func- 
tion b(x, k s = 1). Note, that the wave- numbers q, param- 
eterizing the Fourier transform of tp, i> = J dx^(q)e lqx , 
are linked with quasi-momenta k as q — k ± 2n. Here 
and below n = 0, 1, 2, . . . numbers the Brillouin zones. 
For values of k s ^ 1 we have observed the quasi-solitonic 
pulses, which, while travelling, leave behind the trail of 
small amplitude radiation, Fig. 2(a). The radiation ef- 
fect becomes noticeably stronger for k s — > k . Spectra of 
the radiating solitons have a distinct peak, see Fig. 2(b), 
which is absent in the spectra of the ideal resting solitons. 
The overall results of the extensive series of numerical 
experiments unambiguously indicate that solitary pulses 
moving through the optical lattice continuously emit ra- 
diation with certain spectrally localized quasi-momenta. 
By analogy with electromagnetic case we term this radi- 
ation as transition radiation of matter waves. 

The initial conditions used above in the form of 
the sech envelope superimposed on the Bloch function 
g s {x)e %ksX are difficult, though probably not impossible 
to prepare in the real experiments. However, current ex- 
perimental techniques allow straightforward preparation 
of the gaussian matter-wave packets with spectrum cen- 
tered around the zero quasi-momentum and setting lat- 
tices in motion. The lattice moving with velocity 2k s will 
then effectively shift the central momentum of the wave 
packet to k s . 2k s equals to the group velocity of the 
free, i.e. without the lattice, matter waves with quasi- 
momentum k s . In turn, the backward scattering of mat- 
ter waves from the moving lattice is expected to create 
the second strong and other peaks in the spectrum of the 
wave-packet. The results of modelling of Eq. (1) with 
moving potential cos2(x — 2k s t) and initial conditions 
in the form of the gaussian packet are shown in Figs. 
2(c,d). Taking into account shift of the axes, one can see 
that these results are in remarkable agreement with those 
obtained using the solitonic initial conditions, see Figs. 
2(a,b). Good resolution of the radiation peaks in Figs. 
2(b,d) suggests even that they can possibly be recorded 
experimentally. 

To understand and give analytical interpretation of the 
observed radiation we develop a perturbative approach, 
allowing us to predict quasi-momenta and the amplitude 
of the emitted wave and, thereby allowing an estimate of 
the decay rate of a solitary pulse. To proceed we form 
the ansatz 

V> = i> s {x,€,t) + ip(t,x), (5) 

where the first term approximates the solitonic part of 
the wave function, ip s — g s (x)A s (^,t)e lksX , and the sec- 
ond one is an arbitrary perturbation. Substituting (5) 
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into (1) and assuming that \ip\ is small we find that evo- 
lution of ip is governed by 



idt<p + L<p-2\il) a \ 2 <p-il%<p* = S(x,t), 



(6) 



where S(x,t) = —idtips — Lxp s + \ip s \ 2 ip s is the source 
term, which is different from zero because ip s is not an 
exact solution of Eq. (1). 

The structure of the radiation tail observed in the nu- 
merical modelling corresponds to the spatially extended 
lattice eigenmode. The natural mechanism for excitation 
of a selected eigenmode from the continuum is the energy 
and wave-number resonance with one Fourier component 
of the source term S. The resonance condition ensures 
that there exists a lattice mode which is always in-phasc 
and therefore interferes constructively with one Fourier 
component of the source term. The latter can be rep- 
resented in the form S = e iksX ~ iest ~ iKt Y,j hj{x)fj(^), 
where the sum is taken over all the terms appearing 
in the right-hand side after substitution of the explicit 
expression for the soliton and hj(x) are some func- 
tions with a period of ir. Replacing /j(£) through their 
Fourier integrals fj(£) — J dQfj(Q)e l ®t and recalling 
that £ = x — e' s t, one can easily find that the lat- 
tice modes satisfying the resonance condition are those, 
which have quasi-momenta k = k r obeying the condition 
e(k r ) = k + e s + e' s Q r . Q r — Q(k r ) is the detuning of 
the wave-number of the resonant wave q r = fe r ± 2n from 
the soliton quasi-momentum k s , i.e. Q r = k r ± 2n — k s . 
Using that e(k r ) = e(q r ) the resonance condition can be 
rewritten in the form 



e(q r ) = k + e s + (q r - k s )e' s . 



(7) 



The geometrical meaning of Eq. (7) is clear. The right- 
hand side of (7) equals to the energy of the Fourier com- 
ponent of the soliton (4) , while its left hand side is simply 
the energy of the linear dispersive wave. Eq. (7) can be 
solved for q r by plotting the tangent line to the peri- 
odically extended dispersion characteristic at the point 
e ~ e s . Then one should make the parallel upshift of this 
tangent by k and find points of the intersection with the 
dispersion characteristic itself. Topologically it is self- 
evident that only resting solitons do not produce any 
real roots of Eq. (7) and do not radiate into the lattice 
modes. For e' s ^ Eq. (7) has infinitely many real roots 
corresponding to different values of n. Root q r for n = 0, 
i.e. when q r = k r , as a function of k s is shown in Fig. 
1(b). Dots on this graph show (k r ,k s ) pairs measured 
from the direct numerical modelling of Eq. (1). Mod- 
elling of Eq. (1) has not revealed any resonances in the 
Brillouin zones with n / 0. It indicates that coupling 
into the higher-order resonances is negligible, primarily 
because spectral strength of the source term is very weak 
for the values of Q r with n ^ 0. Note, that the secondary 
less intense spectral peaks, seen in Figs. 2(b,d), are de- 
scribed by the peaks in the Fourier spectra of b(x, fc r ,s), 



and are given by q TtS — fc r>s ± 2n. At the same time, 
multiple roots of Eq. (7) are not linked by any simple 
algebraic expression. 

In order to calculate the amplitude of the radiated 
wave we form the anzats 



ip ~ ip s (x,£,t) +g r (x)e 



(8) 



where subscript V refers to the quantities calculated at 
k = k r . Substituting (8) into Eq. (1) we take into ac- 
count that Lg r (x)e ik <- x W(t,£) ~ -g r (x)e ik ' x £ r W(£,t). 
The latter expression can be easily inferred by compari- 
son with Eq. (3). Assuming that Sk s <C \k s — k r \ <C 2, 
where Sk s is the width of the soliton in the quasi- 
momentum space and 2 is the width of the Brillouin zone, 
and disregarding terms nonlinear in W, one gets the av- 
eraged equation for the amplitude of the radiation field: 
s = id t W + C r W, where s= (P, -P*) T , W = {W,W*) T , 
P = e ms-kr)[a i {£ s -£ 2s }R s+ a i \R s \ 2 R s ], D r [-id ( ] = 
—e r + ie' s di= + £ r [—idt\, 



-D r + a 2 \R s \ 2 onR 
-aXR 2 e 2l ^ k -- k ^ Dt 



2 p 2i£(k 3 -k r ) 

*M 2 



(9) 



ai = / dxg 2 s gf, a 2 = 2 / dxg 2 \g r \ 2 , a 3 = a s J dxg s g* - 
J dx\g s \ 2 g s g;, and a 4 = J dxg s g*. 

The approximate solution we obtain for W is 

w(t, ~ -iC{e(C0 - e(a + (t(e' s - <))} . (io) 

Here £ = sgn(e' r — e' s ), is the Heaviside function and C 
is the amplitude, which can not be generally expressed 
in a closed analytical form and, therefore, was calculated 
numerically. C characterizes the spectral intensity of the 
source term for k — k r . Heaviside functions in Eq. (10) 
describe the tail of the radiation field having the length 
t\e' r — e' s \. The tail starts at the soliton (£ = 0) and 
extends beyond (e' r < e' s ) or in front (e' r > e' s ) of it. 

Substituting the ansatz (8) into the conservation law 
dt f dx\ip\ 2 = 0, we can estimate the rate, 7, of the 
transfer of particles from the soliton to the radiation: 
7 = C 2 \e' s — e' r \. Plot of log 10 7 vs k s for two values of 
k is shown in Fig. 3(a). Naturally, the rate of transfer 
increases, when detuning \k r — k s \ decreases for k s — > fco, 
see Fig. 1(b). This is because for k s — > fc the lattice 
mode is resonant with the most intense central part of 
the soliton spectrum. Contrary, \k r — k s \ increases for 
k s — ► 1, and the radiation amplitude decays almost ex- 
ponentially. For example, semi-analytical Eq. (10) gives 
that the initial condition with k s = 0.87 used to generate 
Fig. 2 has the initial decay rate 3 x 10 4 particles per 
second. Providing that this rate is constant in time, the 
soliton life-time would be ~ 0.033s. However, radiation 
carries away both density and momentum from the soli- 
tonic part of the field. Therefore solitonic parameters n 
and k s are in fact functions of time. In particular, we 
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have found that the radiation results in convergence of 
the soliton quasi-momentum to some limit value, which 
is always closer to 1 than the initial k s . Thus, radiation 
emission slows the soliton. For example, taking initial 
conditions with n = 0.01 and k s = 0.87 we have observed 
that over 0.05s (1000 dimensionless time units) k s shifts 
to ~ 0.887, see Fig. 3(b), and k to ~ 0.006. The rela- 
tively small increase of k s is accompanied by the decrease 
of the soliton decay rate, which drops, accordingly with 
plots shown in Fig. 3(a), from 3 x 10 4 to ~ 10 2 particles 
per second. Let us stress again that the soliton life-time 
tends to infinity when k s — > 1. Correspondingly, the es- 
timated soliton life-time increases dramatically to ~ 10s. 
Fig. 3(b) also shows the numerically computed slowdown 
in the growth of the number of particles in the radiation 
component of the condensate. Soliton decay rate inferred 
from these data and corresponding theoretical points cal- 
culated for instant values of the soliton parameters k and 
k a are in good agreement, see inset in Fig. 3(b), which 
confirms validity of our theoretical method. 

In summary, we have reported transition radiation by 
matter-wave solitons moving through the optical lattice. 
This effect extends family of the already known and re- 
lated quantum radiative effects such as, e.g., sound emis- 
sion by precessing quantized vortices [16] and by dark 
matter- wave solitons oscillating in a harmonic trap [17]. 
Note, also that our main conclusions and techniques can 
be used to predict and analyze radiation by spatial op- 
tical solitons moving in nonlinear photonic crystals [12]. 
This work was partially supported by the INTAS project 
211-855. 
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FIG. 1. (a) Energy- momentum diagram for the linear mat- 
ter-waves in optical lattice: n — Brillouin zone is shown, (b) 
Dependence of the radiation quasi-momentum, k r , from the 
soliton quasi-momentum, k 3 . Full diagonal line in (b) corre- 
sponds to k r = k s . Dots in (b) mark (k s ,k r ) pairs measured 
from the modelling of Eq. (1). Dashed vertical lines in (a) 
and (b) mark the ±fco points with e" = 0. j3 = 0.33 




4 



FIG. 2. Results of the direct numerical modelling of Eq. 
(1). (a,c) Square roots of the atomic density as functions of 
x. The inset in (a) shows fine details of the spatial profile of 
the radiation. (b,d) Corresponding Fourier spectra. (a,b) are 
obtained for the initial solitonic wave packet with k s = 0.87 
and k — 0.01. (c,d) are obtained for the initial gaussian wave 
packets with the zero central momentum and lattice moving 
with velocity 2k s — 1.74. Integration time t = 700 corre- 
sponds to 0.035s. The physical number of particles in the 
soliton shown in Fig. 2(a) is ~ 10 3 , its velocity is3xl0~ 3 m/s 
and width at the half-height is ~ 4/im. 
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FIG. 3. (a) Logarithm of the particle transfer rate from the 
solitonic to the radiation part of the wave function as a func- 
tion of the soliton quasi-momentum for k = 0.01 and 0.006. 
Dashed vertical line marks the e" = point, (b) Temporal 
evolution of the normalized number of particles N ra d in the 
radiation component of the field (left axis) and corresponding 
dynamics of the soliton quasi-momentum (right axis). N ra d is 
calculated as integral of \ip\ 2 over the tail behind the soliton. 
Estimate for physical number of particles in the radiation tail 
is given by 3 x 10 3 N ra d- Inset shows the soliton decay rate 7 
as function of t: solid line - numerical modelling and dots - 
theoretical results. 
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